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ABSTRACT 

With the imminent promise of constraints on the epoch of reionization from observations 
of microwave background anisotropies, the question of whether or not the standard Cold 
Dark Matter (CDM) model permits early reionization has been subjected to detailed 
investigation by various authors, with the conclusion that reionization may occur at 
quite high redshift. However, it is widely accepted that this model is excluded, as when 
normalised to the COBE observations it possesses excessive galaxy clustering on scales 
below tens of megaparsecs. We examine the trends of observations, first in a fairly 
model independent way, and second by considering variants on the standard CDM 
model introduced to resolve the observational conflicts. We conclude that the epoch of 
reionization favoured by the observational data is typically considerably later than the 
standard CDM model suggests, and amongst models which may fit the observational 
data only the introduction of a cosmological constant leads to a reionization redshift 
close to that of standard CDM. 
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1 INTRODUCTION 



O O ne of the most important tests that a model of large scale 
I structure must pass is the Gunn-Peterson test, which indi- 

^ cates that the universe is nearly completely reionized by a 
* redshift of five or greater (Gunn & Peterson 1965; Schneider 
et al. 1989). Until recently, there has been little prospect 
of imposing an upper limit on the ionization redshift, but 
that is now promised via measurements of microwave back- 
ground anisotropies on degree scales and less; if reioniza- 
tion is sufficiently early then the expected 'Doppler' peak 
on scales of around a degree may be sharply suppressed due 
to the photons scattering off the reionized electrons at lower 
redshifts (Sugiyama et al. 1993). The ever increasing cata- 
logue of anisotropy experiments already hints strongly that 
the peak is present at around the correct height for mod- 
els based on adiabatic perturbations without reionization, 
arguing against very early reionization (Scott et al. 1994). 

Reionization provides a fascinating link between differ- 
ent regions of the power spectrum of density perturbations 
(Couchman & Rees 1986), since it is normally assumed that 
it is induced by early formation of massive stars which pro- 
duce an ionizing flux. The first objects to form are usu- 
ally assumed to be galaxies of mass around 1O 6 M0, corre- 
sponding to the mass enclosed in a comoving region of di- 
ameter around 0.01 Mpc, but their effect on the microwave 
background alters the interpretation of anisotropies which 
originate from perturbations on scales of hundreds of mega- 



parsecs (Vittorio & Silk 1984; Bond & Efstathiou 1984). 
Reionization will also raise the Jeans' mass which may in- 
fluence the formation of structures somewhat larger than 
the ones which induce the reionization (Couchman & Rees 
1986). 

Recently, two detailed papers have appeared, one by 
Tegmark, Silk & Blanchard (1994) (hereafter referred to 
as TSB) and one by Fukugita & Kawasaki (1994) (here- 
after FK), which examine the question of reionization in 
the standard Cold Dark Matter (CDM) model of structure 
formation. The standard CDM model comprises a spatially 
flat universe dominated by cold dark matter, containing 
baryons with density consistent with standard nucleosyn- 
thesis (though many authors assume zero baryon density), 
and with an adiabatic, gaussian, scale-invariant spectrum 
of density perturbations. We shall define precisely what we 
take as the standard CDM model later. [TSB, and in ad- 
dition a follow-up paper by Tegmark & Silk (1994), also 
investigated several other models. However, they did not 
normalise their spectra to the COBE observations; the final 
result is extremely sensitive to the normalisation.] Although 
the approaches of TSB and FK were rather different, the for- 
mer concentrating on estimation and the latter on numerical 
simulation, very similar results were obtained. The conclu- 
sion was that with reasonable parameter values, the CDM 
model will reionize early, at a redshift around 35; optimistic 
parameter choices obviously lead to yet earlier reionization. 
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This is easily early enough to have an important impact on 
degree scale anisotropies (Sugiyama et al. 1993). 

However, it is widely accepted that standard CDM is ex- 
cluded by large scale structure observations, because when 
normalised to COBE it possesses excessive power on scales 
below some tens of megaparsecs. Since the epoch of reioniza- 
tion is very sensitive to the amount of power in the spectrum 
on very short scales, the suspicion is therefore that observa- 
tions favour substantially later reionization. It is our aim in 
this paper to take advantage of the detailed work by TSB 
and FK to study the implications of large scale structure ob- 
servations for the reionization epoch, at first making as little 
reference to specific models of structure formation as we can, 
and then discussing features specific to the different ways in 
which one can conceive of extending standard CDM to sat- 
isfy the observations. In cold dark matter dominated mod- 
els with critical density, the analysis is particularly simple; 
other models feature more complicated evolution of density 
perturbations at moderate redshift which must be taken into 
account. Important aspects of the results we present are that 
throughout we shall normalise the spectra to COBE, using 
the recent results of Bunn, Scott & White (1994), and that 
we shall incorporate the effect of baryons into the transfer 
functions we use. 



2 THEORY AND OBSERVATIONS 

2.1 Standard CDM and its variants 

All of the models considered here assume that large scale 
structure originates as an adiabatic gaussian density per- 
turbation, whose spectrum at horizon entry is more or less 
scale independent with any scale dependence parametrised 
as a power of the comoving wavenumber k. In order to 

have any chance of agreeing with observations these models 
require non-baryonic cold dark matter. 

In general, such models are specified by the spectral 
index n, the Hubble parameter H whose present value is 
parametrised as 100/ikms -1 Mpc -1 , the nature and density 
of the non-baryonic dark matter, and the value of the cosmo- 
logical constant. The density of baryons as a fraction of the 
critical density, Qb, is determined by nucleosynthesis to be 
n B h 2 = 0.013±0.002. The standard CDM model has n = 1, 
h = 0.5, pure cold dark matter, total density O ma t = 1 
and no cosmological constant. Normalised to the COBE ob- 
servation it gives too much power on smaller scales, so a 
number of variants of it are currently under consideration. 
These are conveniently characterised by modifying just one 
assumption of the standard CDM model, though it is im- 
portant to recognise that nature may have chosen to modify 
at least two of them. The options under active investigation 
at present are 

(i) The tilted* model (Bond 1992; Liddle et al. 1992; Cen 
et al. 1992) which has spectral index n < 1. 

* In an inflationary context, the introduction of tilt is accom- 
panied by the possibility of a gravitational wave contribution to 
COBE (Davis et al. 1992b; Liddle & Lyth 1992). This will always 
make reionization more recent, and we shall neglect this possibil- 
ity here. 



(ii) The low Hubble constant model (Bartlett et al. 1994) 
which has h < 0.5. 

(iii) The Mixed Dark Matter (MDM) model (Bonometto 
& Valdarnini 1984; Shan & Stecker 1984; Davis et al. 
1992a), which replaces some of the CDM by hot dark 
matter (normally assumed to be a massive neutrino) 
with density fi„. 

(iv) The cosmological constant model which replaces some 
of the CDM by a cosmological constant corresponding 
to a vacuum energy density V ac- 

(v) The open universe model which simply throws away 
some of the CDM. 

2.2 Strategy for calculating the reionization 
redshift 

Following TSB, we estimate the reionization redshift in three 
stages. 

First we need the fraction / of baryons that has to be 
bound into stars in order to provide sufficient ionizing flux 
that the reionized fraction is unity. This is uncertain by or- 
ders of magnitude. However, the joy of examining reioniza- 
tion in the context of models based on gaussian random 
fields is that this fraction is an exponentially sensitive func- 
tion of the amplitude of perturbations and hence of redshift, 
resulting in an uncertainty of only a modest factor in calcu- 
lating the latter. 

Throughout we shall use the estimates of / derived by 
TSB. They write / as a product of several factors each 
of which is estimated separately, and conclude that the 
'Middle-of-the-Road' estimate for the collapsed fraction re- 
quired to induce complete reionization is/~8xl0~ 3 . The 
optimistic and pessimistic estimates (where in their termi- 
nology 'optimistic' refers to earlier reionization) are given as 
/ ~ 4 x 10~ 6 and / ~ 0.8 respectively. Note that the pes- 
simistic estimate is very pessimistic indeed, being a product 
of pessimistic estimates for four separate factors (likewise for 
the optimistic estimate). The largest individual uncertainty 
is the fraction of hydrogen burned within a short time of 
star formation. Once reionization has occurred, the ioniz- 
ing flux which continues to be emitted is enough to prevent 
subsequent recombination (Fukugita & Kawasaki 1994). 

The next step is to suppose that this fraction / is given 
by /(> M m i n , Zi on ), where /(> M, z) is the fraction ^ of the 
mass of the universe that is bound into objects with mass 
bigger than M at redshift z, and M m i n is the minimum 
galaxy mass at the redshift of reionization. Two criteria are 
of relevance for estimating this minimum mass. The bary- 
onic Jeans' mass governs pressure support against collapse, 
and the mass must be large enough that the cooling time, 
which governs dissipation, is no greater than the dynamical 
time (Blanchard et al. 1992). These two criteria are broadly 
similar, and the normal assumption is that M m i n is around 
1O 6 M0, with leeway of around an order of magnitude in 
either direction. 

The final step is to estimate /(> M,z) within a given 

t In the hope of avoiding confusion, we have eschewed the usual 
notation Q(> M, z), as Q conventionally gives the fraction of a 
critical density, whereas for reionization we are interested in the 
fraction of the total matter density, regardless of its value. 
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model using the Press-Schechter approximation (Press & 
Schechter 1974). This states that /(> M, z) is equal to twice 
the fraction of space in which the density contrast, smoothed 
on the scale M, exceeds a threshold 6 C . Since the probability 
distribution of the filtered density contrast is gaussian, this 



gives 



/(> M, z) = erfc 



V2a(M,z) 



(1) 



where a(M,z) is the dispersion of the density contrast 
smoothed on a scale M at redshift z, and 'erfc' is the com- 
plementary error function. 

To complete the definition of the Press-Schechter ap- 
proximation one has to specify the filter function used to 
perform the smoothing, and the threshold 6 C .^. We shall use 
a top-hat filter, and take the threshold value 6 C = 1.7 which 
is motivated by a spherical collapse model. A gaussian filter 
is also commonly employed. For a given M the gaussian- 
smoothed dispersion is higher than the top-hat-smoothed 
dispersion by a significant scale-dependent factor (Liddle & 
Lyth 1993). This means that the equivalent value of 6 C is 
significantly lower, the equivalence being of course only ap- 
proximate because the relative factor is scale-dependent. 

The Press-Schechter approximation is expected to be 
reasonable provided that the mass fraction, or equivalently 
the dispersion, is significantly less than unity. In other words, 
it is expected to be valid in the linear regime, where the fil- 
tered density contrast is evolving according to linear cosmo- 
logical perturbation theory. TV-body simulations can check 
its validity (as long as the mass fraction is not too small), 
and also provide an empirical 'best fit' value for 6 C . Results 
vary (Efstathiou & Rees 1988; Brainerd & Villumsen 1992; 
White et al. 1993; Lacey & Cole 1994; Ma & Bertschinger 
1994), but such simulations typically suggest that the ap- 
proximation is at least roughly correct, with S c lying within 
twenty percent or so of the theoretically motivated 1.7. As 
one expects, the values of S c suggested for Gaussian smooth- 
ing are typically somewhat lower than those suggested for 
top hat smoothing. 

With pure CDM and critical density, the evolution of 
(t(R, z) in the linear regime is simply 



a(R,z) = a(R,0)/(l + z). 



(2) 



where a(R,0) is the linearly evolved quantity at the present 
epoch. Substituting this expression into (1) gives 



1 + Zion 



^a(M min ,0) 



erfc _1 (/) . 



(3) 



This is the formula we use for the optimistic and 
Middle-of-the-Road estimates of /. With the pessimistic es- 
timate / ~ 0.8, most of the mass is supposed to have col- 
lapsed before reionization occurs. In that case the quasi- 
linear theory underlying the Press-Schechter approximation 
ceases to be reliable, and we replace it by the crude es- 



Neither TSB nor FK are very specific about the type of smooth- 
ing. TSB use a top-hat filter but quote some lower values of S c 
which are at least in part lower because the papers cited used a 
gaussian filter. FK use a gaussian filter but take S c — 1.69 which 
is normally motivated by top-hat collapse. 



timate that a significant fraction of mass collapses when 
a(M m i n , z) ~ 1, leading to 

1 + Zion = a(M min ,0). (4) 

2.3 The transfer function 

The top-hat filter of comoving radius R is given by 
w/i D^ siakR cos kR 

w{kR) = TkW~TkW' (5) 

where k is the comoving wavenumber. The relation between 
mass and radius for the top-hat filter is 

3 



M = 1.16 x 10 12 fi mat h' 1 



R 



(6) 



Following the notation of Liddle & Lyth (1993), the present 
dispersion is 



a(R,0) = 9H { 



Jo 



k 4 T 2 (k)S 2 H (k)W 2 (kR)^. (7) 
k 



The transfer function T(k) is normalised to unity as k — » 0. 
The initial spectrum Sj 1 (k) is defined here to be independent 
of k for a flat (n = 1) spectrum; for a tilted spectrum one 
has &%(k) oc A;" -1 . The usual quantity P(k) is proportional 
to k6%(k). 

The transfer function for a given version of the CDM 
model can be calculated by numerical integration of the rel- 
evant evolution equations. Unfortunately, published calcu- 
lations are not completely adequate for our purpose and we 
shall have to do a certain amount of improvisation. For the 
moment let us consider the case where all the dark matter 
is cold. 

Both TSB and FK take the transfer function from 
Bardeen et al. (1986) (henceforth BBKS) without further 
comment. However, there are two important questions. 
Firstly, calculation of transfer functions to very small scales 
is troublesome because of the difficult integration over many 
oscillations of the baryon fluid; most published transfer func- 
tions are only accurate down tolfe 1 Mpc at best ( see the 
discussion in Liddle & Lyth (1993)). Although the BBKS 
transfer function has the correct asymptotic form including 
the logarithmic correction to the large k behaviour, they 
supply no indication of the accuracy of its amplitude on 
scales as short as we are interested in. Secondly, the BBKS 
transfer function is for zero baryon density (for which the 
usual Qh scaling is exact); however, baryons will reduce the 
short scale power, and even at 8/i _1 Mpc this effect is known 
to be around 10% for h = 0.5. 

There is little we can do about the former problem, but 
we can at least attempt to allow for the latter via an em- 
pirical scaling of the transfer function with baryon density 
as publicised by Peacock & Dodds (1994). On intermedi- 
ate scales this brings the BBKS transfer function into good 
agreement with calculations including baryons (at least for 
the range of Qb we are interested in); that the BBKS trans- 
fer function then possesses the right asymptotic form, un- 
like other parametrisations, may to some extent address the 
question of accuracy at short scales too. We therefore take 
the transfer function as 

In (1 + 2.34g) 



Tcdm(<?) 



2.34<j 
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[l + 3.89g + (W.lg) 2 + (5.46g) 3 + (6.71g) 4 ] ~ 1/4 , (8) 
where q = fc/[fi mat /i 2 exp( — 2Qb) Mpc -1 ]. 

2.4 Observations 

The above estimate of the reionization redshift needs as in- 
put the value of (t(R, 0) on the scale R ~ 0.01 Mpc. Observa- 
tions which have been used to constrain large scale structure 
theories do not reach down to such small scales, but they do 
cover four orders of magnitude from roughly the size of the 
observable universe down to about one megaparsec, with 
fairly tight constraints available on various scales. Most of 
the data can be summarised as a value of (t(R, 0) which al- 
lows it to be presented in a single plot, Figure 1. A full dis- 
cussion of the data requires a separate paper (Liddle et al. 
1994) [see also Lyth & Liddle (1994)], and here we shall only 
summarise it with particular emphasis on the crucial COBE 
normalisation. Our summary shall concentrate on O ma t = 1; 
a full generalisation to arbitrary O ma t will be given elsewhere 
though we shall mention some aspects later. 

Over the desired range, the dispersion varies by several 
orders of magnitude while we are interested in differences 
of a factor two. Consequently, we shall plot all the data as 
its ratio to the prediction of standard CDM; since standard 
CDM certainly fits all the data to within a factor two or so 
(though not within the actual errors), this greatly clarifies 
the graphical presentation. 

In Figure 1 the error bars indicate something like a 1-<t 
range, but the upper and lower limits indicate rather firm 
results. 

(i) The COBE measurement of microwave aniso- 
tropics. The most accurate measurements of the mi- 
crowave background anisotropy presently available are 
those of the COBE satellite. The analysis and inter- 
pretation of these measurements remains an evolving 
field. A very elegant analysis was recently carried out 
by Gorski et al. (1994). Using the Sachs- Wolfe ap- 
proximation for the C; (the expected value of the l- 
th anisotropy multipole), they found that a good fit 
is obtained with a spectral index n anywhere in the 
range 0.6 to 1.4, and that all fits gave essentially the 
same value for Cg. For n = 1 they concluded that 
the expected quadrupole <3 rms -ps = (19.9 ± 1.6)jttK. 
Very recently Bunn et al. (1994) have done an inde- 
pendent fit of the Sachs- Wolfe approximation, finding 
the value Q rms -ps = (21.1 ± 1.6)jttK which is slightly 
but not significantly higher. However, they go on to 
fit to a full calculation of the predicted anisotropics, 
including the tail of the Doppler Peak. More or less 
independently of the nature of the dark matter and 
the values of h and Qb, they find 

Q rms _ PS (n) = (19.9 ± 1.5) exp [0.69(1 - n)] fiK , (9) 

and this is the normalisation we shall adopt. The re- 
lation between &h and <3 rms -ps can be derived ana- 
lytically for power law spectra (Liddle & Lyth 1993) 
and for n = 1 is 6 H = 2.26 x 10 _5 (Q rmB _p S /19.9/iK). 
We shall schematically represent the COBE point as 
being at 4000/i~ 1 Mpc, though when discussing spe- 



cific models they shall be normalised using the Bunn 
et al. result directly, 
(ii) The galaxy correlation function. This is very nat- 
urally given in terms of a(R,0). Data are now suffi- 
ciently copious to provide constraints over more than 
an order of magnitude in R. We take the values given 
by the compilation of data sets by Peacock & Dodds 
(1994), converted to a(R,0) according to their pre- 
scription. In Figure 1 they are represented by a band, 
whose edges correspond to 1-sigma errors excluding 
overall normalisation. The larger error bars at each 
end of the band indicate the uncertainty in the over- 
all normalisation, which allows one to shift the entire 
data set up or down while preserving its shape. The 
data is scaled for the best fit bias parameter with O ma t 
held at unity, bi = 0.85^2^, - This is determined from 
redshift distortions and nonlinear effects in the spec- 
trum (Peacock & Dodds 1994). However, the following 
data more stringently constrain the overall normalisa- 
tion. 

(iii) Peculiar velocity flows. The bulk flow around us 
smoothed on a scale 40/i~ Mpc is quite well measured 
(Bertschinger et al. 1990). The velocities are actually 
sensitive to longer scales in the dispersion (t(R, 0), and 
in practice this measurement samples a(R,0) on a 
scale of around 90k' 1 Mpc (Lyth & Liddle 1994). Un- 
fortunately, as it is a single measurement the 'cosmic 
variance' is high, dominating the observational errors. 
An analysis using a comparison of peculiar veloc- 
ity flows with IRAS galaxy distributions (Dekel et 
al. 1993) is more useful, giving bi = Q.T^q® at 
95% confidence level. Using the upper limit to nor- 
malise the Peacock and Dodds data gives the limit 
a(8/i _1 Mpc,0) > 0.6 sh own in Figure 1. 

(iv) Galaxy cluster abundance. This directly measures 
the dispersion on a scale around 8h~ Mpc, using 
either the Press-Schechter approximation or numer- 
ical simulations. From the calculations of White, Ef- 
stathiou & Frenk (1993), we deduce that if the viri- 
alised mass of Abell clusters with richness class > 1 is 
less than 1.5/T 1 x 1O 16 M , then cr(8/i _1 Mpc, 0) < 1.0. 
This mass is three times as big as the one estimated 
from temperature and virial velocity measurements, 
though we note that it might be consistent with grav- 
itational lensing measurements. In any case, it seems 
reasonable to regard the above bound as rather firm, 
and it is shown in Figure 1. 

(v) Early object formation. The spectrum must con- 
tain enough power on short scales to reproduce the ob- 
servations of quasar abundances (Efstathiou & Rees 
1988) and the fractional density in damped Lyman 
alpha systems (Mo & Miralde-Escude 1994; Kauff- 
mann & Chariot 1994). This provides lower limits on 
a(R,z), where z is the relevant redshift; results can 
then be scaled to z = 0. The lower limits we show 
assume purely cold dark matter but are not much dif- 
ferent in viable MDM models. They should be rather 
firm as they correspond to very loose interpretations 
of the measurements (Lyth & Liddle 1994; Liddle et 
al. 1994). 

For our final table of results, it is useful to charac- 
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Figure 1. Observational values of the present linearly evolved dispersion of the smoothed density contrast. These values apply to critical 
density CDM and MDM models, but not to low density models. They are plotted as a fraction of the predicted value for the standard 
CDM model. From the right: the solid triangle represents COBE; the solid square is from the bulk velocity flow; the band indicates the 
1-sigma allowed region of the galaxy correlation points according to Peacock & Dodds (1994) (15 points are plotted and the jaggedness 
is noise) while the error bars at each end of the band indicate the uncertainty in bias from this method; the upper and lower limits at 
8/i -1 Mpc, coming from cluster abundance and bulk flows respectively, are limits on the normalisation at high confidence which are much 
stricter than those from the galaxy correlation function alone; the other lower limits are from the quasar abundance (right) and from 
the amount of gas in damped Lyman alpha systems (left). Finally, the notches on the left indicate reionization redshifts for the standard 
CDM model (top value), and modifications of it which invoke tilt and/or a low Hubble constant to reduce the small scale power. The 
estimates of z; on are positioned at the scale corresponding to M — 10 6 Mq, because the reionization redshift depends only on a(R, 0) on 
that scale. The aim is to extrapolate through the observations to the notches from which the reionization redshift can be directly read 
off. For MDM models the observational data are almost unchanged, but the reionization redshifts must be rescaled as indicated in the 
text. For open and cosmological constant models this Figure does not apply at all; they must be treated separately. 



terise the fit to the data by two parameters, which indicate 
how well COBE normalised models fit the galaxy correlation 
function. One is the overall normalisation, which we char- 
acterise by a(8/i -1 Mpc, 0). As discussed its allowed range 
is 

a(8/i _1 Mpc,0) = 0.6 to 1.0, (10) 

if O ma t = 1. The other parameter we need is some measure 
of the shape. In models without a hot component, we follow 
Peacock & Dodds (1994) in defining the shape parameter by 

T = fi/iexp(-2fi B ) - 0.32 - lj . (11) 

In MDM models we estimate an effective T from the slope of 
the dispersion at 25/i -1 Mpc. At the 2-sigma level the range 
of r that fits the data of Peacock & Dodds (1994) is 

T = 0.22 to 0.29 , (12) 

(see also Kofman et al. (1993)). 



3 THE EPOCH OF REIONIZATION 

In this section we calculate the reionization redshift, first in 
the standard CDM model, which is ruled out by observation, 
and then in variants of it which may be consistent with 
observation. The results shall be displayed in Table 1. 

3.1 Standard CDM 

In the standard CDM model the reionization redshift is 
given by (3) and (4). Using the COBE normalisation, one ob- 
tains a(M min ,0) = 33.9, 29.1 and 24.6 for M min = 1O 6 M , 
1O 6 M and 10 7 M Q respectively. If we take the Middle-of- 
the-Road estimate for the collapse fraction / required to 
reionize, we get z- lon = 53, 46 and 38 for the three values 
of M m i n . On the other hand, keeping M m i n fixed, the opti- 
mistic and pessimistic estimates for / multiply these figures 
by 1.5 and 0.6 respectively. As this latter is the larger un- 
certainty (and also dominates the uncertainty in 6 C ), from 
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now on we shall fix M m i n to be 10 M@. Accounting for our 
different approach regarding the normalisation of the spec- 
trum and the incorporation of baryons, these results agree 
with TSB, unsurprisingly as we have followed their estima- 
tions. The principle reason for the somewhat higher redshift 
is the higher normalisation to COBE. 

From Figure 1 it is clear that standard CDM is ruled out 
by observation, so we go on to consider variants of it, retain- 
ing always the COBE normalisation. In each case we want to 
choose the additional free parameter to give a reasonable fit 
to the slope and magnitude of the galaxy correlation data. 

3.2 Tilted and low Hubble constant models 

For the 'tilted' and 'low Hubble constant' variants of the 
standard CDM model, perturbation growth is just as in stan- 
dard CDM, and so (3) and (4) remain valid. Consequently, 

1 + ^ion is just reduced in proportional to a(M m i n , 0). The 
aim in altering these parameters is to try and fit the large 
scale structure observations; extrapolating the data shown in 
Figure 1 through to the notches, which have been placed at 
the scale corresponding to M m i n = 10 Mq, gives the reion- 
ization redshift. Without worrying about particular param- 
eter values, with the strong downward trend enforced by the 
galaxy correlation data it seems hard to imagine that these 
variants can give a reionization epoch to be much above 

2 = 20 in these variants, even with these data points raised 
collectively upwards as much as allowed. 

In Table 1 we give the normalisation and shape parame- 
ter for the representative choices h ~ 0.3 and n ~ 0.7. Either 
of them does a fair job of satisfying the data and they both 
give reionization around z- lon = 10. 

3.3 Mixed Dark Matter models 

In an MDM model, a fraction fi„ of the matter is in the 
form of hot dark matter. Its free-streaming leads to a scale 
dependent suppression in the growth of perturbations, but 
for the short scales relevant to reionization the effect is sim- 
ply a slowed growth rate across all relevant epochs, with the 
hot component playing no part in the gravitational insta- 
bility (Davis et al. 1992a). The growth rate during matter 
domination is 

vmdm(M,z) oc (l + z)- a , (13) 
with 

a = (V25 - 24fi„ - l) /4 . (14) 

Consequently, if we know the value of o"(10 6 Mq,0), it is 
simple to track it back to determine the reionization redshift. 

Unfortunately, there do not exist transfer functions for 
MDM which are accurate, or even have the right asymptotic 
form, for the scales we are interested in. Instead, we follow 
Davis et al. (1992a) and use a sudden transition approxima- 
tion. This assumes that in an MDM model the perturbation 
growth up to the end of the radiation era can be taken as the 
same as CDM; the transition to matter domination is then 
taken as instantaneous and the suppressed growth law above 
is assumed to take effect immediately. If accurate transfer 
functions are obtained in the future, it will be interesting to 
compare them with this approximation. 



A subtlety to be taken into account is that in the stan- 
dard incarnation of the MDM model, where the hot com- 
ponent is a massive tau neutrino, this neutrino becomes 
non-relativistic before matter domination, which then oc- 
curs earlier than in the usual case where the tau neutrino 
acts as radiation. The usual epoch z eq ~ 24000Q ma t/i 2 is 
replaced by z eq ~ 27000fi ma t/i 2 , since each neutrino species 
contributes 12% of the present radiation density. 

We shall quote results for two values of the hot com- 
ponent density. The originally favoured fi„ = 0.3 (giving 
a = 0.8) has lost favour due to its apparent inability to ex- 
plain the amount of gas in damped Lyman alpha systems 
(Mo & Miralde-Escude 1994; Kauffmann & Chariot 1994; 
Ma & Bertschinger 1994), and so we shall also quote for 
Q„ = 0.15 (giving a = 0.91). In Table 1 we give the val- 
ues of a(8/i -1 Mpc, 0) and T for these two choices, calculat- 
ing the former from the transfer functions of Schaefer and 
Shafi (1994) (which are accurate on these scales), and esti- 
mating an equivalent value for the latter from the slope of 
0"mdm(-R, 0) at R = 25/i _1 Mpc. Taking into account the 8% 
uncertainty in the COBE normalisation, both choices give 
acceptable values for a(8/i -1 Mpc, 0), and the values of T 
are also reasonable though an intermediate choice would do 
better. 

Most of the data in Figure 1 applies at redshift zero, 
and the correction to the object abundances is small for 
reasonable fi„. Consequently, Figure 1 still indicates the 
trends in the data. However, the reionization values given 
by the notches must be rescaled, with 1 + z- lon replaced by 
(1 + Zion) 1 /". Alternatively, the approximation above gives 
an analytic formula for the central reionization redshift 

l + z lon = 47 1/ "(27000/i 2 ) 1_1/ ". (15) 

For fi„ = 0.3, the central redshift estimate is z- lon = 13, 
while for fi„ = 0.15 it is Z[ OIl = 27. The growth suppres- 
sion is playing some role in making these values larger than 
they would otherwise have been, though particularly in the 
case of the latter reionization is earlier through the actual 
normalisation of the power spectrum being rather high. The 
slower perturbation growth also results in the spread in red- 
shifts from the uncertainty in / being larger than in other 
cases. 

3.4 Cosmological constant and open universes 

Finally we come to the cosmological constant and open mod- 
els, where the present matter density O ma t is less than one. 
In these models the COBE normalisation is amended and 
all observational data need reassessed, so Figure 1 cannot 
be used at all. We content ourselves by examining specific 
parameter values, which we choose to be O ma t — 0.3 and 
h ~ 0.8. These do as well as any in minimizing the disagree- 
ment with observation (but note that for the open case the 
present age of the universe would be only 10 Gyrs, while 
with a cosmological constant it is 12 Gyrs). 

Before a redshift of at most a few, O ma t is close to unity 
and the growth of perturbations is as in the critical density 
case. As a result the predicted reionization redshift is the 
same for a given choice of the spectrum 8%{k) of the pri- 
mordial density perturbation. Afterwards the growth slows 
down, so that by the present the density contrast is sup- 
pressed by a factor g relative to the critical density case. 
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Table 1. Estimates of the reionization epoch Z[ on , assuming a minimum collapse mass of 1O 6 M0. The uncertainty in this mass is 
dominated by the uncertainty in the collapsed baryon fraction required to induce reionization. For information we include the predicted 
values of the parameters o"(8/r — 1 Mpc, 0) and T, which measure the normalisation and slope of the galaxy correlation function. From 
observation the target range of V is 0.22 to 0.29. For critical density models the target range of cr(8/i- 1 Mpc, 0) is 0.6 to 1.0, but for the 
low density models in the last block it is 0.9 to 1.4. 



Model 


h 


n 


^mat 




a(8/i _1 Mpc,0) T 

low estimate 


Zion 

best guess 


high estimate 


Standard CDM 


0.5 


1.0 


1.0 





1.41 


0.45 


28 


46 


69 


Tilted CDM, n = 0.7 


0.5 


0.7 


1.0 





0.71 


0.31 


6 


10 


16 


CDM, h = 0.3 


0.3 


1.0 


1.0 





0.64 


0.23 


8 


13 


20 


MDM, fi„ = 0.15 


0.5 


1.0 


1.0 


0.15 


1.08 


0.27 


16 


27 


44 


MDM, fi„ = 0.30 


0.5 


1.0 


1.0 


0.3 


0.93 


0.19 


7 


13 


22 


CDM, fi vac = 0.7 


0.8 


1.0 


0.3 





1.64 


0.23 


25 


39 


61 


Open CDM 


0.8 


1.0 


0.3 





0.48 


0.23 


11 


19 


29 



For the matter density contrast the suppression factor is 
accurately parametrised by (Carroll et al. 1992) 

g = f Q mat [0l 7 t - fi vac + (l + ^f) (l + ^)] "!(16) 

To use this result we need the COBE normalisation. A com- 
plete discussion along the lines of Bunn et al. (1994) has not 
been given, but the literature contains sufficient information 
that we can extract a reliable normalisation. 

For the cosmological constant model there is no spa- 
tial curvature, so the Fourier expansion can still be used 
to define the spectrum 6 H (k). As in the standard CDM 
model we take it to be scale independent, though note that 
the most detailed comparison of this model with observa- 
tions (Kofman et al. 1993) included some tilt. With this 
assumption, the COBE data have been fitted by Bunn & 
Sugiyama (1994); they find that the overall fit to COBE is 
quite poor (see also Sugiyama & Silk (1994)), and quote a 
best fit Qrms-ps = 21.6/ttK. Unfortunately they do not give 
the corresponding value of Sh, but we have calculated it 
by making the reasonable assumption that the Sachs- Wolfe 
approximation is good for the quadrupole. Because the po- 
tential fluctuations are induced by Sp rather than Sp/p, the 
matter fluctuations are then bigger by a factor l/fi ma t. The 
correction due to the time dependent potential induced by 
the cosmological constant has been calculated by Kofman & 
Starobinsky (1985); for a given Sh, the expected quadrupole 
is about 9% bigger than in the critical case for the value 
fimat = 0.3. This increase happens to coincide with the 
increase in the best fit value of <3 rms -ps, so we conclude 
that the COBE normalisation of Sh is larger just by the 
factor 1/0.3. The central reionization redshift is therefore 
Zion = 39, the closest to the standard CDM value of all the 
models we have studied. Including the growth suppression 
factor g = 0.78 gives the other results in Table 1. 

For the open universe model, there is no unique gener- 
alisation of the concept of scale-invariance in the spectrum 
on scales comparable to the curvature scale. This makes the 
COBE normalisation ambiguous. To break the ambiguity, 
we take the view that the scale dependence should be that 
given by the inflationary prediction using the conformal vac- 
uum (Lyth & Stewart 1990). 



Using the inflationary shape for the spectrum, Sugi- 
yama & Silk (1994) have calculated the C;. For O ma t = 0.3 
they find that the shape mimics the critical density Sachs- 
Wolfe prediction with a spectral index n ~ 1.4. As men- 
tioned earlier, Gorski et al. (1994) found that the COBE 
data can be fitted with 0.6 < n < 1.4, so this value is 
marginally allowed, though as in the cosmological constant 
case the fit will be poor compared with the critical density 
case. We therefore normalise the open model by requiring 
that the predicted value of Cg agrees with the value found 
in the fit of Gorski et al. To determine this normalisation, we 
use values of C; given by Kamionkowski et al. (1994) (they 
actually normalise to the ten degree variance of the temper- 
ature anisotropy, but they also give the Ci). The resulting 
normalisation of Sh is 0.49 times that for critical density, 
leading to a later reionization redshift z- lon ~ 19. Further, 
the suppression factor of g = 0.45 is stronger than for the 
cosmological constant model, leading to a much lower value 
for a(8/i _1 Mpc,0). 

Concerning whether these models fit the large scale 
structure data, we note that the observations must be rein- 
terpreted. The shape is still given by T from (11), but the 
observational value of o"(8/i _1 Mpc, 0) increases. Peacock & 
Dodds (1994) give the scaling of their power spectrum with 
O, and White et al. (1993) give the scaling of the cluster 
abundance. In combination, these give a target range 

a^^T 1 Mpc, 0) = 0.9 to 1.4 (17) 

in either open or cosmological constant models with O ma t = 
0.3. We shall make a more detailed comparison of these mod- 
els with observation elsewhere. 



4 CONCLUSIONS 

The most important aspects of the results that we have pre- 
sented here are firstly that we have normalised all of the 
models to COBE, and secondly that we have taken account 
of baryonic corrections to the transfer functions. The nor- 
malisation to COBE is at the high end of the possibilities 
considered by TSB for standard CDM, but is considerably 
lower than that they considered for variant models. The in- 
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elusion of baryons reduces the result for z- lon by about 20% 
in typical models. For O ma t = 1, we have presented a col- 
lection of observational data across a wide range of scales 
in a suitable form for extrapolation towards the scales rele- 
vant for reionization. It is clear that observations favour an 
epoch of reionization considerably later than that found in 
the standard CDM model. 

We then examined a series of COBE normalised models. 
Most options favoured by the data give reionization around 
a redshift ten to twenty. However, the cosmological con- 
stant model, with its significantly higher normalisation of 
fluctuations, gave an answer considerably higher than the 
rest and almost as high as the standard CDM model. The 
other type of model capable of giving rather early reioniza- 
tion is an MDM model with a low choice for the density of 
the hot component; this arises via a combination of mod- 
est growth suppression and the model's natural tendency to 
explain the shape of the galaxy correlation function with a 
higher a(8/i -1 Mpc, 0) than other models. 

The most important hurdle that models must surpass 
in connection with reionization is the Gunn-Peterson test 
(Gunn & Peterson 1965). For some variant CDM models, 
this could conceivably be tricky if the most pessimistic sce- 
nario is realised, but the question clearly cannot be ad- 
dressed with existing understanding. Lower limits on the 
short scale perturbation spectrum from early object for- 
mation such as damped Lyman alpha systems (Mo & 
Miralde-Escude 1994; Kauffmann & Chariot 1994; Ma & 
Bertschinger 1994) promise to be much more secure. 

The other role that reionization can play is in the damp- 
ing of degree scale microwave background anisotropies (Vit- 
torio & Silk 1984; Bond & Efstathiou 1984). This damping 
has been analysed recently by Sugiyama et al. (1993), and 
proves to be very sensitive to the reionization redshift. At 
z- lon ~ 50, towards the top end of the standard CDM predic- 
tions, the Doppler peak may be damped almost completely 
away. Even if Z[ OIl is around 20, which appears more feasible, 
they find a suppression of over 30% which is certainly signif- 
icant and seems likely in the cosmological constant model. 
By a redshift of 10, appropriate to several models, the effect 
is very small. 
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